Tree structure of the percolating Universe 
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We present a numerical study of topological descriptors of 
initially Gaussian and scale-free density perturbations evolv- 
ing via gravitational instability in an expanding universe. We 
carefully evaluate and avoid numerical contamination in mak- 
ing accurate measurements on simulated fields on a grid in a 
finite box. Independent of extent of non linearity, the mea- 
sured Euler number of the excursion set at the percolation 
threshold, 5c, is positive and nearly equal to the number of iso- 
lated components, suggesting that these structures are trees. 
Our study of critical point counts reconciles the clumpy ap- 
pearance of the density field at 5c with measured filamentary 
local curvature. In the Gaussian limit, we measure \5c\ > a 
in contrast to widely held belief that |5c| ~ o", where cr^ is the 
variance of the density field. 



The large scale structure of the universe can be char- 
acterized topologically at various scales by studying the 
excursion sets of the density contrast (5(x); i.e., the sets 
= {x|<5(x) > <5™} or E^^ = {x|^(x) < <5™}. An 
interesting choice of (5th in either cases is the percola- 
tion threshold, ^c, for which the excursion set includes 
an infinite, connected structure ||^. 

Accurate determination of be is a challenge both from 
the theoretical point of view - no analytic calculation of 
(5c exist so far in the three-dimensional case; the exper- 
imental point of view - numerous spurious effects, such 
as discreteness, finiteness of the sampled volume, etc., 
can contaminate the measurements, and the results de- 
pend on the method employed Here, we proceed 
as in Ref. [Q and consider site percolation on a grid, i.e., 
neighboring sites above/below the density threshold are 
connected through faces. 

We analyze the topology of Ef^^ in terms of the three 
dimensional local curvature of (5(x) given by the Hessian 
matrix, d^S/dxidxj. The number of negative eigenval- 
ues, /, of the Hessian matrix at any point classifies the 
local density structure into one of four distinct classes: 
region with / = 3 is in a "clump" (c); / = 2 in a "fil- 
ament" (f); / = 1 in a "pancake" (p) and / = in 
a "void"(v). In particular, this classification applies to 
critical points of the field, where the gradient dbjdxi = 0. 
Connectivity of the excursion set E'f^^ is dictated solely 
by the number, C/, of critical points of each class / within 
Ef^^. For clarity, we replace the numeral index / of C/ 
by the alphabetic label / = 0, 1, 2, 3 ^ v, p, f, c. 



At a simple qualitative level, the role of critical points 
in outlining the connectivity of excursion sets can be un- 
derstood quite intuitively. Connectivity at percolation 
is along special field lines ("ridges" and "river beds") 
threading the critical points. Filament saddle points (f- 
saddles) are thus crucial for percolation in overdense re- 
gions {E'^^J. Indeed, a large fraction of them lie along 
ridges connecting two local maxima. When Sth > <5 at 
a saddle point, the point drops out of the excursion set, 
thus disconnecting the corresponding clumps. Analogous 
reasoning can be applied to percolation in underdense 
regions where a large fraction of pancake saddle points 
(p-saddles) lie on river beds that connect local minima. 

In detail, the situation can be more complicated than 
this simple vision, e.g., saddle-points connected through 
ridges or river beds. However, the work of Morse rigor- 
ously establishes a simple link between the number distri- 
bution of critical points (local property related to curva- 
ture) and global connectivity (topology) of the excursion 
sets of any smooth field lending credence to our above 
intuitive line of argument A simple, but powerful, 
relation exists between the Euler characteristic, X @, of 
Ef_^^ and a linear combination of the four critical point 
counts, Cj. It can indeed be shown that despite all the 
possible complexity Q, 

X± = ± [ Ce - Cf + Cp - Cv] . (1) 

The observed structures in the galaxy distribution are 
believed to have evolved through gravitational instability 
from small initial fluctuations, usually postulated to be 
a Gaussian random field. We pay particular attention to 
the Gaussian limit since at large scales the memory of 
initial conditions is retained in the evolved system. At 
small scales, non linear clustering creates an asymme- 
try between overdense and underdense regions possibly 
changing the properties of the percolating network and 
the value of Sc- To study clustering in the non linear 
regime we use A^-body simulations. 

In this paper we have two goals: (i) obtaining an accu- 
rate measurement of 6c as a function of smoothing scale; 
(ii) characterizing the topology of the percolating net- 
work through a combination of the Euler number, critical 
points and local curvature. We attempt to exploit eq. (||) 
to understand the global properties, such as percolation, 
of the evolving density contrast in the universe in terms 
of the distribution of (local) critical points. 
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In our numerical experiment, we consider iV-body sim- 
ulations (described in |^]) with scale- free Gaussian initial 
conditions, where the power-spectrum as a function of 
wavenumber, k, spectral index, n, and amplitude. A, is 
given by P{k) — Ak", n = —2,-1,0. This covers the 
observationally relevant range of values in the scaling 
regime 1 ^ ^ ^ 100 Mpc. The simulations were generated 
with a particle-mesh code |^ with 256'^ particles in a peri- 
odic cubic box of size L. A mesh resolution N,- = 256 was 
used to compute the forces. The density field was com- 
puted using Cloud-In-Cell interpolation Q on grids of 
various sizes A^g = 64, 128, 256 and then convolved with 
Gaussian windows of scales £/L = 1/Ng, 2A/Vg, . . . , 0.04, 
0.08 to insure sufficient differentiability [Q. The corre- 
sponding variance of the smoothed density distribution 
will be denoted by a'^{i). To increase the dynamic range 
of our measurements, we take advantage of self-similar 
behavior for scale-free initial conditions ||ll[] and combine 
several time snapshots of each simulation with the appro- 
priate scaling {L,£) {L/s,£/s), where s is the linear 
theory correlation length {cTfi^gg^,.{s) = 1). For Gaussian 
smoothing, s = [a^A r[(n 3)/2]/(47r2)]i/(»+3)^ where a 
is the expansion factor of the Universe, arbitrarily fixed 
to unity at the beginning of the simulation. 

To study the Gaussian limit, a set of 10 pure Gaussian 
realizations of the density field, (5(x), at each grid size and 
spectral index, {Ng,n), were generated and smoothed at 
the same set of scales, (/L, as in the A^-body samples. 

The gradient and the Hessian of the density field are 
computed by least square fitting a second order hyper- 
surface to the values 6 computed at each site and its 26 
closest neighbors. 

We have explored the wide range of smoothing scales 
£/L and box sizes Ag to thoroughly quantify the numeri- 
cal effects contaminating topological properties. The op- 
timal set of parameters (A^g, £/L) = (256, 0.01) is adopted 
to minimize the following contamination effects (partly 
discussed in [Q): 

(i) The grid must be sufficiently thin compared to the 
smoothing scale, Ng£/L > 2.5. With this constraint, 
our procedure for computing numerical derivatives of 6 
provides extremely accurate results. 

(ii) Periodic boundaries enhance percolation in a finite 
box but the effect is negligible if £/L ^ 1. On the other 
hand, finite box size suppresses large scale [C'(i)] power, 
which can be viewed as an increased effective spectral 
index. From our Gaussian samples, we find that combi- 
nation of these two competing effects is inconsequential 
if £/L < 0.01. 

(iii) Inherent limitations of the A^-body approach: spa- 
tial resolution of the forces imposes N^£/ L > 1. To evade 
effects of A^-body relaxation and unrealistic non linear 
couplings between the large 0{L) modes, it is required 
that 2/Afg/^ <s/L and s/L <{).!, respectively. Finally, 
discreteness effects are insignificant since an 8-fold dilu- 



tion of the number of particles does not change the results 
substantially. 

First column of panels in Fig. |] shows the measured 
percolation threshold, (5c, as a function of scale, £/ s. In 
agreement with previous works [|l2|,^, the asymmetry 
due to gravitational clustering between overdense (-I-) 
and underdense regions (— ), J^t > \5'^\^ increases with 
the level of non linearity and with —n. It is very small 
for n = 0. For underdense regions, the asymptotic regime 
(5(7 ^ —1 is naturally reached when £/s <C 1. The value 
of 5'^ I (T is maximum near transition to the nonlinear 
regime (see also Ref. |^]), £/s ^ 1, where initial struc- 
tures are optimally enhanced by large-scale motions , 
and is larger than the Gaussian limit, the difference in- 
creasing with —n. As noted later, 5'^ /a can become small 
in the highly nonlinear regime due to lost coherence. 

In the Gaussian limit, we find slight increase in 5c with 
-n; \5c/(t\ = 1.6±0.1, 1.4±0.1, 1.3±0.1 for n = -2,-1,0, 
respectively. Error bars reflect the dispersion in the 10 
realizations and possible contaminations (i) and (ii) dis- 
cussed above. These results should be contrasted with 
that of Ref. ||l^, where \5c/<j\ ~ 1 is obtained indepen- 
dently of n. However, the measurements in Ref. JT^ use 
a smaller grid, Ag = 64, and might be affected by ef- 
fects (ii) mentioned above. The other difference is that 
the field is softened by a power-spectrum cut-off at small 
scales instead of our Gaussian smoothing. 

Second column of panels in Fig. ^ displays the fraction 
// of space occupied by clump, filament, pancake and 
void regions in Kf^ . As shown earlier from a qualitative 
point of view pq |, the simulated universes at percola- 
tion are seen to be composed mostly of filamentary and 
clumpy regions, in the respective fractions of 55% and 
35 — 40% in the Gaussian limit. The filamentary nature 
becomes dominant in the non linear regime and the re- 
sults do not depend significantly on n. 

Third column of panels in Fig. Q shows the measured 
fraction of critical points of each kind. Note that maxima 
are dominant, explaining the rather clumpy appearance 
of the largely filamentary (in terms of curvature) perco- 
lating network in Fig. |^. In fact, dumpiness augments 
with non linearity as the fraction of maxima increases 
and explains the rather small 5'^ /cr despite increased fil- 
amentarity when tr ^ 1. 

Fourth column of panels in Fig. ^ shows the "reduced" 
Euler number 8 = x/Ctot, Ctot = X^/C/, as a function 
of density threshold measured at various stages of the 
simulations, corresponding to a = 1, 32, 128 for n = 0, — 1 
and a = 1, 16, 32 for n = —2. The value of a increases 
with the abscissa of the left local maximum of 8. For 
a = 1 , the measurements agree almost perfectly with the 
analytic prediction in the Gaussian limit as expected. 
The 8 curves retain overall shape with time, except that 
they are progressively compressed near the origin and 
stretched at large ^th/c. This perhaps indicates that 
non linear clustering does not significantly influence the 
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FIG. 1. Topological properties of scale-free universes with Gaussian initial conditions, (a) The measured percolation threshold 
Sc as a function of scale. Each panel corresponds to the spectral index n indicated on it, respectively n = —2, —1, from top to 
bottom. Open (filled) circles correspond to measurements in the A'^-body simulations for overdense (underdense) regions. The 
errorbars represent deviations from self- similarity. The rightmost square gives the Gaussian limit obtained from the average over 
10 random realizations. The solid curve is l/a(£/s), corresponding to the expected asymptotic regime 5^ — > —1 for percolation 
of underdense regions in the limit £/s 0. (b) The measured fractions, //, of clump/filament /pancake/void regions in E+ 
as functions of scale (respectively, open and filled circles, stars and asterisks). Again, the rightmost points give the Gaussian 
limit, where the //'s can also be computed analytically as functions of Sth ||15|. (c) As in (b) but the fractions of the four types 
of critical points, Ci/Ctot- (d) The measured Euler characteristic as a function of density threshold at three different stages of 
the simulations, as explained in the text (solid curves). The number of isolated clusters Ndus/Ctot is shown by dots. The solid 
curve that corresponds to initial conditions matches perfectly the analytic Gaussian limit Q (dashes) . Symbols with horizontal 
error bar indicate the measured value of the percolation threshold, 5c. 



topology of the system. 

There are three regimes dictated by the sign of the first 
derivative of £. With increasing (Jth/o", (i) at small Sjn/a, 
mostly local minima drop out of E"^^ , creating cavities, 
thus increasing £; (ii) at intermediate Sth/(J, p-saddles 
start to drop out too, and cavities connect together, thus 
the value of £ decreases; (iii) in the large Sth/ct region 
even f-saddles drop out, breaking up the ridges to create 
isolated clusters, thus £ increases again. At this last 
stage, the contribution of minima and p-saddles to the 
Euler number can be largely neglected, x ~ Cc — Cf . 

Fourth column of Fig. 1 embodies the main result of 
the paper. It shows that percolation occurs at positive 
values of the Euler characteristic x > in regime (iii), 
and specifically at x — -^cius, where A'cius is the num- 
ber of "isolated clusters" (disconnected regions of Eg'^^). 
This means that at percolation threshold, isolated clus- 



ters, including the largest percolating one, have close to 
tree structure [0, i.e. they contain no or few loops. In- 
deed, in this case each cluster contributes unity to the to- 
tal Euler number adding up to x ~ .^cius- The f-saddles, 
the only type available in number in the regime (iii) can 
produce only loops, which decrease Euler number of indi- 
vidual clusters, and being abundant would decrease the 
total X much below A'cius, which is not what is observed. 
Thus we conclude that the percolating network is a tree. 
Fig. H gives a visual presentation of a percolating tree. 

In summary, we used very simple but novel statistics in 
cosmology - critical points of the density field and its local 
curvature to characterize the topological properties of the 
universe at percolation. The most important results of 
our analyses are the following: 

(a) At percolation, the Euler characteristic is positive, 
equal to the number of isolated clusters suggesting that 
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FIG. 2. 3-D view of the overdense percolating network in 
the n = — 1 simulation, at a = 128 and smoothed at £ = 0.02. 
It consists of spherical-like high density regions dominated 
by 'clump' points (/ — 3) connected by bridges built mainly 
from points with filamentary local curvature (7 = 2). No 
significant number of closed overdense loops is present. 

most of these have tree structure (with almost no loops). 
This result is valid for overdense as well as undcrdcnse 
excursions. In particular, we find in the Gaussian limit 
that If^d/cr > 1, at variance with common belief (e.g. 
Ref. iQ) that site percolation in Ef^^ should occur at 
|5c|/c = 1, where the Euler number of E^^ is zero. 

(b) In terms of curvature, the overdense structures at 
percolation are mostly filamentary and clumpy in agree- 
ment with previous works Filamentarity in- 
creases with level of non linearity. However, maxima are 
the dominant critical points in the excursion set. With 
isotropic smoothing, this explains its rather clumpy ap- 
pearance as illustrated by Fig. ^. 

Our efforts to identify and address sources of contam- 
ination in our numerical experiments suggest that it is 
actually difficult to measure the percolation threshold ac- 
curately in a real galaxy catalog, where one of the main 
issues is contamination by discreteness effects. These 
impose a large smoothing scale, a fact difficult to rec- 
oncile with the requirement i/L ^ 0.01, even for current 
or forthcoming large redshift surveys such as the Two- 
degree Field Q or the Sloan Digital Sky Survey [p^ . 

However, it is still possible to measure local curvature 
and critical point counts, because these statistics should 
not require as small values of £/ L to be measured with 
sufficient accuracy. Although, here we focus on critical 
point counts at the percolation threshold, these counts as 
functions of 5th are of interest. In particular, the statis- 
tics Cj and // can be computed explicitly as functions 
of (5th in the Gaussian limit [^. As the genus of the 
boundary of Ef^^ [Q, they could be used in real galaxy 
catalogs to constrain the nature of the primordial den- 
sity field. To do that properly, several issues remain to 
be addressed, such as the anisotropics induced by redshift 



distorsion in three-dimensional galaxy catalogs or the ef- 
fects of biasing between the observed galaxy distribution 
and the underlying total matter distribution. 
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